Elastic lines on splayed columnar defects studied numerically 
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We investigate by exact optimization method properties of two- and three-dimensional systems of 
elastic lines in presence of splayed columnar disorder. The ground state of many lines is separable 
both in 2d and 3d leading to a random walk -like roughening in 2d and ballistic behavior in 3d. 
Furthermore, we find that in the case of pure splayed columnar disorder in contrast to point disorder 
there is no entanglement transition in 3d. Entanglement can be triggered by perturbing the pure 
splay system with point defects. 
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I. INTRODUCTION 

Transport properties of type-II superconductors are in- 
fluenced by the presence of various kinds of disorder . 
Pinning of vortex lines hinders their motion, which, in 
response to an apphed current, causes dissipation. From 
the practical point of view it is highly desirable to avoid 
the appearance of vortex creep which gives rise to a fi- 
nite resistivity. It was proposed by Hwa et al. that 
splayed columnar defects resulting from heavy ion irra- 
diation of superconducting samples, would significantly 
enhance the vortex pinning, and thus reduce the vortex 
creep leading to a higher critical current density jc. An 
additional important aspect might be the inhibition of 
vortex motion due to forced entanglement induced by 
the disorder 

The predictions concerning jc have been been veri- 
fied in experiments on samples with different sources of 
splayed defects Kwok et al. @ reported that 

well above the matching field B^, where the density of 
vortices equals to the density of columnar defects, the 
irreversibility line of sample with splay is below of the 
one without splay defects. As an evidence of strong en- 
tanglement above the matching field it was discovered 6] 
that although the irreversibility is decreased in samples 
with splay defects values of jc are still greatly increased 
compared to unirradiated samples. On the other hand 
the comparison of samples with one and two families of 
splay defects did not reveal any differences in the values 
of jc 0- Molecular dynamics simulations for B > 
show the increase of jc in the samples with splay disor- 
der H, 0- The authors of 9] have also performed simu- 
lations on samples with B < B^ for which they obtained 
a milder enhancement of jc- On the basis of these ob- 
servations they suggested that in the case B > B^ the 
additional increase of jc is due to vortex entanglement. 

The ground state properties of an ensemble of flux lines 
in such disordered environments has never been investi- 
gated before. Single flux line properties in the presence 
of tilted columnar defects at zero temperature were stud- 
ied by Lidmar et al. They show that the behavior 
of the lines depends on the energy distribution of the 
lines. This is manifested in roughening, or mean-square 



displacement as a function of sample height i?. 
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\H) = {v{zf)^-{v{z))l^H^^. 



(1) 



Here z is the distance along the height direction, and 
r the transverse displacement. According to Lidmar et 
al. the roughness exponent Q is sensitive to the shape of 
the distribution of the tilt angle and the energy distribu- 
tions of the defects For instance with an opening 
angle of 90° and for a uniform energy distribution the 
roughness exponent in 2d is ^ = 3/4, in contrast to the 
point-disorder result 2/3 In samples with a fixed 
starting point a single line has the following geometry. 
It occupies a splay defect until a jump to a energetically 
more favorable one takes place. The lines undergo jumps 
from splay defect to splay defect so that the average dis- 
tance between two successive jumps grows as Az ~ z [lol| . 
Thus though the jump density decreases with growing z 
the roughening exhibits a non-trivial scaling. 

The natural question arises how the single-line physics 
outlined above changes in the presence of many inter- 
acting lines, at a constant line density p. In this pa- 
per we study this in both two (2d) and three dimensions 
(3d). There are studies of the role of disorder in two- 
dimensional samples [l^ . while the three-dimensional 
case corresponds to bulk superconductors. We find that 
at finite line densities the physics changes in particular 
as the roughening of lines is concerned: the roughness 
exponent becomes C = 1/2 in 2d and ^ = 1 in 3d. The 
ballistic behavior of the lines leads to the absence of true 
entanglement in 3d. 

This paper is organized as follows. In Sec.^we explain 
the details of our numerical model. The scaling of rough- 
ness is discussed in Sec lIIII We also study in Sec. IIVI if 
inserting a new line causes rearrangements in the con- 
figuration of previously inserted lines, ie. as the density 
is increased. The results shown in Sec. demonstrate 
that in pure splay disorder lines do not entangle. Entan- 
glement can be induced by perturbing splay with point 
disorder. The summary and discussion are presented in 
Section EH 
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FIG. 1: Schematic figure of the graph consisting of splayed 
columns in two and three dimensions. In the latter case note 
the two additional bonds connecting the columns. 



II. MODEL 

We consider the following model of a system of in- 
teracting lines in a two- or three-dimensional disordered 
environment: The lines live on the bonds of a graph con- 
sisting of an ensemble of splayed columns embedded in 
a box with a width L and a height H (Fig. QJ. Each 
column is described with a transverse coordinate r(z) at 
height z from the bottom level: 



r(z) = rg + az 



(2) 



where rg is a random point on the basal plane and a is 
a randomly chosen variable that defines the magnitude 
and the orientation of the tilt of a given column. 

The number of columnar defects M is set to M — L 
and M = in two and three dimensions respectively. 
The average distance between two nearby columnar de- 
fects defines the in-plane length unit of the model. The 
lines enter the system at the bottom plane (z = 0) and 
they exit at the top plane (z — H), and wc use both fixed 
and free entrance points, and free exit points. The case 
of a fixed entry point and a free exit point is that consid- 
ered by Lidmar et al. for a single line lOj, who discuss an 
experimental scenario for the same. Also rough enough 
surface could pin the end points of flux lines. Strong 
pinning along large surface steps has been reported in 
Ref. [13. More recent experimental results [lj| indicate 
that surface pinning have measurable effects on the flux 
lattice dynamics. 

In this paper we study the case where the columns 
have a random, uniformly distributed orientation within 
a cone with a fixed opening angle. In our model changing 
the opening angle of the cone is equivalent to a rescaling 
of the system height, for which reason we consider only 
the opening angle of 90°. 

The number N of lines threading the sample is fixed 
by the prescribed density p — N/M. Within this version 
of the model number of lines cannot exceed the number 
of columns, which means in the case of splay disorder 
that B < B^. The graph could be in principle modified 
such that lines can traverse the system also between the 
columnar defects which would correspond to _B > B^. 



In the transverse directions we use both periodic and 
open boundary conditions. In the latter case the defects 
crossing the boundary arc cut such that lines cannot fol- 
low them across the system boundaries. Note that one 
can not let flux lines escape from truly open boundaries 
since the line density would decrease with z, the lon- 
gitudinal coordinate. In order to avoid correlations in 
systems with H > L we first construct a graph of size 
iJ'' from which we cut a piece of size L. 

In 2d lines can change defects only at their crossings 
and in 3d when the defects are close enough to each other. 
In 3d columns are connected by introducing an extra seg- 
ment between the lines at the shortest mutual distance 
whenever this distance is shorter than a fixed value Tj,. 
When Tc is kept relatively small compared to the aver- 
age distance between the columnar defects the choice of 
the energy cost of the bond connecting two defects can 
be arbitrary. We choose zero energy cost and Tc = 1.4. 
We have also checked the value rc — 0.5 which did not 
change the results. 

We model the disordered environment by assigning a 
random (potential) energy Cij = r^M to each bond (i, j), 
where Vij is the Euclidean length of the bond and u is 
a random number which determines the type of disor- 
der. In the case of splay disorder the random variable u 
is drawn independently for each column such that bonds 
along a given column share the same value of u. Point dis- 
order is modeled with a similarly constructed graph. The 
only difference is that u is drawn independently for each 
bond which destroys the correlations along the columns. 
For point disorder we use uniformly distributed u and for 
splay disorder the following distribution in order to make 
comparisons to single line results by Lidmar et al. ^10] : 



Piu) 



(3) 



for < w < 1 and otherwise P{u) — 0. With = 1 this 
reduces simply to a uniform distribution. 

For computational convenience we focus on the short 
screening length limit (15(] . We restrict ourselves to hard- 
core interactions between the lines, which means that 
their configuration is specified by bond-variables Xij G 
{0, 1}, Xij = 1 indicating that a line segment occupies a 
bond between nodes i and j, and Xij = indicating that 
no line segment occupies this bond. The total energy of 
the line configuration is given by 



(4) 



where the summation is performed over all bonds. The 
corresponding Hamiltonian in a continuum limit is given 
by the following formula: 



n 



dr. 



dz 



+ Vr[v.,{z),z] 



+ V^nt[r.(z)-r,(z)]} , 



(5) 
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FIG. 2: Schematic figure of adding lines to the graph, (a) 
The ground state of a single line with arrows indicating the 
direction of the optimal path, (b) The bold line is the optimal 
path from bottom to top in the residual graph. The residual 
graph is constructed by changing the energy costs of the bonds 
to negative along the dotted line. In addition, the direction 
of paths in the residual graph on the dotted bonds must be 
opposite to the direction of the previously inserted path shown 
in (a), (c) The final line configuration after the path segments 
with opposite directions are canceled. 



where ri(z) denotes the transversal coordinate at longi- 
tudinal height z of the i-th flux hne. The interactions 
T4nt[ri(z) — are hard core repulsive and the dis- 

order potential yr[ri(z),z] is (5-correlated in the case of 
point disorder and strongly correlated along columns in 
the case of columnar disorder. The elastic energy is mim- 
icked in our numerical model by the positivity of all en- 
ergy costs per unit length, which ensures for instance the 
minimization of the line length in the absence of random- 
ness. 

At low temperatures the line configurations will be 
dominated by the disorder and thermal fluctuations are 
negligible. Therefore we restrict ourselves to zero tem- 
perature and focus on the ground state of the Hamilto- 
nian Q. Computing the ground state now corresponds 
to finding TV non-overlapping directed paths traversing 
the graph along the bonds from bottom to top. One 
has to minimize the total energy of the whole set of the 
paths and not of each path individually (n.b.: already 
the two-line problem is actually non-separable For 
a single line problem the lowest energy path can be 
found straightforwardly with Dijkstra's shortest path al- 
gorithm. In the case of many lines Dijkstra's algorithm 
is applied successively on a residual graph instead of the 
original graph as illustrated in Fig.|21 In the residual 
graph the properties of occupied bonds are modified as 
follows: If Xij = 1 we set e-°^"*^'^' = —eij, i.e. negative, 
and require that the direction of a new path must be 
in the opposite direction with respect to the direction of 
previously inserted paths. The path segments with two 
opposite directions on a given bond are canceled. For 
details see [l7| . 

Although the lines cannot occupy the same bond of 
the lattice they may touch in isolated points if the nodes 
in these points have more than three neighbors. This is 
now the case only in 2d as exemplified in Fig^ Since 
we want to calculate the roughness of lines we need to 
determine the individual lines, for which we use a local 
rule. In 2d the line identification is unambiguous if we 



simply require that the lines cannot cross. 

In our model flux lines are confined to the defects with 
no possibility to enter the bulk (except in 3d if the jumps 
from one defect to another are counted as such) . We also 
tested the case where a homogenous bulk between the 
splayed defects is represented by a set of densely packed 
columnar defects with zero tilt and with a constant en- 
ergy cost per unit length. Here, after a flux line has found 
a splayed defect with a lower energy cost per unit length 
in z-direction compared to the one in the bulk one recov- 
ers the pure splay disorder behavior. Thus, including the 
possibility for lines to travel also in the bulk introduces 
another crossover length. We have checked numerically 
that this scenario holds. After a cross-over system size, 
which depends on the ratio of the energy costs in the bulk 
and on the defects, for single line systems we obtain the 
roughness exponents of pure splay disorder. 

In many line systems with a flxed ratio between the 
number of splay defects and the number of inserted flux 
lines it is not guaranteed that all lines can find a de- 
fect which is energetically more favorable than the bulk. 
Depending on the energy cost of the bulk there can be 
only a fraction of splayed defects which are energetically 
more favorable compared to the bulk. If the number 
of such defects is smaller than the number of fiux lines 
there will be a fraction of lines fbuik that stay in the bulk 
throughout the whole sample. As a function of the en- 
ergy cost per unit length in the bulk Ubuik this scales like 
fbuik = g(ubuik/ P°') with g{x) a scaling function that de- 
creases monotonically from 1 at a; = to zero for x = 1 
and a depending on the shape of the energy distribution 
of the splayed defects (Eq.Q). 

Since the lines staying in the homogenous bulk have 
no transverse fluctuations our central results on the geo- 
metrical properties of flux line systems in splay disorder 
are not expected to change qualitatively. 



III. ROUGHNESS OF THE GROUND STATE 

In this section we focus on the scaling of the average 
roughness w, i.e. the amount of transverse wandering of 
the lines. The mean square displacement the z-th line 
is w'j = {vi{zY)z — {vi{z))1 where Vi{z) is the transverse 
coordinate of the line i at the distance z from the bottom 
level and (r(z))z = jj Jq r{z) dz denotes the average 
along the line from the bottom z = to the top z = H. 
We define as the disorder average of wf averaged over 
all lines. 

In the case of point disorder our model is in the same 
universality class as the directed polymer model accord- 
ing to which the roughness of one line scales as w ~ 
with well known exponents (^2d = 2/3 and Cad « 5/8 in 
two and three dimensions respectively j^J- 2d the 
steric repulsion between the lines leads to collective rear- 
rangement of the lines which yields a logarithmic growth 
of the roughness. In 3d lines can wind around each other 
which suppresses the repulsion resulting in a random 
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walk-like behavior of lines ■ 

For splay disorder, we propose the following simple sce- 
nario. At small z lines do not see each other and exhibit 
single line behavior. Beyond some value of z, which de- 
pends on the density p, the lines cannot further opti- 
mize their configurations and stay on the same defects. 
This leads to a linear growth of the roughness, C = 1 
in 3d due to ballistic behavior. The same arguments on 
the structure of the optimal line configuration hold also 
for 2d. However, in 2d splayed defects can cross each 
other while the individual flux lines are identified such 
that they do not cross the other flux lines. Thus, in 2d 
one has effectively system of hard-core repulsive random 
walks resulting in C = 1/2. Thus one would expect a 
roughness scaling form 



,iL,H)^L'^fiH/L), 



(6) 



where f{x) is a scaling function. In both 2d and 3d, this 
scenario is independent of the splay energy distribution. 

Fig. shows the correctness of this proposition in 3d. 
One can see from Fig. |2Ia) that the roughness of many 
line systems grows linearly w ^ H with no dependence 
on the energy distribution of the defects. According to 
the data collapse in Fig. Ofb) the saturation roughness 
and the saturation height grow linearly with the system 
width L, cf. with Eq.®. 

In 2d there is a collective regime, where the lines ex- 
hibit random walk like -behavior as suggested above. Fig. 
^a) shows that the roughness grows asymptotically like 
w ~ H^^^ for both values of v. The data collapse shown 
in Fig.^b) gives the random walk like scaling also for the 
saturation roughness in agreement with Fig.^Ja). Close 
to the system boundaries the roughness of the lines is sup- 
pressed (due to the way the line identification is made in 
2d) which shows up in the scaling of small system widths 
L. 



IV. SEPARABILITY OF THE GROUND STATE 

We define that the ground state configuration of N 
lines is separable if it can be obtained by adding suc- 
cessively flux lines to the system without modifying the 
configurations of the previous lines. This can be checked 
from the successive shortest path algorithm: whenever 
segments of flux lines are canceled in the residual graph 
this implies that adding a line changed the previous con- 
figuration (see Fig. [SJ. Hence we focus on such seg- 
ments or bonds and calculate the sum of the energy costs 
Es = Cij of such bonds and use it as a measure of the 
separability. When Eg = no flux is canceled due to pos- 
sible rearrangements of line configuration. This means 
that one has fully separable ground state. 

In the case of (splayed) columnar disorder it is obvi- 
ous that with periodic boundary conditions and the full 
freedom of choosing the most favorable starting point 
the lines pick up the defects in the order of their total 
energy cost from the source to the target. Hence, no 
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FIG. 3: Doubly logarithmic plots for the roughness w of 
multi-line systems in 3d. Top: Roughness vs. system height 
H {L = H) for different values of line density p and for two 
defect energy distributions = 1, 2); lines have fixed starting 
points. For comparison are shown also the roughness for point 
disorder (labeled as "p.d."). Bottom: Scaling plot for the 
roughness with fixed system widths for p = 0.2 and v — 1, 
the un-scaled data are shown in the inset. At this line density 
the difference between the data with fixed and free starting 
points is undistinguishable. 



rearrangements of fiux line configurations is needed. In- 
troducing the boundaries or any other distortions to the 
pure columnar disorder reduces the separability of the 
ground state. 

In Fig. 13 we demonstrate the difference between the 
separability of the ground states of splay and point dis- 
order. We consider L x H systems in two dimensions 
with periodic boundary conditions in the direction per- 
pendicular to the 2-axis. In the case of splay disorder 
Eg saturates at a particular system height, because line 
rearrangements do not take place at greater z. Thus Eg 
divided by the system volume goes to zero with increas- 
ing H , implying separability. Compare with the fact that 
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FIG. 4: Double logarithmic plots for the roughness w of 
multi-line systems in 2d. Top: Roughness vs. system height 
H {L = H) for different values of line density p and for two 
defect energy distributions {v = 1, 2); lines have fixed starting 
points. Bottom; Scaling plot of the roughness with fixed 
system widths L and fixed starting points with line density 
p = 0.1 and 2/ = 1. 



Es is linear in H in the case of point disorder, indicating 
non-separability. We calculated also in three dimen- 
sions and observed the same behavior as in the 2d case 
as is depicted in Fig. 



V. ENTANGLEMENT 

In 3d lines can wind around each other resulting in 
topologically non-trivial configurations which we analyze 
here by computing the winding angle of all line pairs as 
indicated in Fig. [7| (c.f. |il9|)- We define two lines to be 
entangled when their winding angle becomes larger than 
2tt . This provides a measure for topological entangle- 
ment 01 , since a topologically entangled pair can not be 
separated by any linear transformation in the basal plane 



FIG. 5; A measure of the separability of the ground state 
in 2d, periodic boundaries, fixed starting points, L=50. Top: 
splay disorder, Es saturates as H grows. Bottom: point dis- 
order, Es grows linearly with H. 

(i.e. the lines almost always would cut each other, if one 
were shifted). 

Sets or bundles of pairwise entangled lines are defined 
such that a line belongs to a bundle if it is entangled at 
least with one other line in the set. In the case of point 
disorder such bundles are spaghetti-like — i.e. topolog- 
ically complicated and knotted sets of one-dimensional 
objects which grow with increasing system height leading 
finally to one giant bundle. In 0] it was shown that for 
point disorder there is a sharp transition from an non- 
entangled phase without giant bundle to an entangled 
phase: The probability Pperc for having an entangled 
bundle of lines that spans the system in the transverse 
direction jumps from to 1 at a critical height He in the 
infinite system size limit (L — > oo), and is described by 
the following finite size scaling form 

Pperc = piL'/^iH ^ H,)) . (7) 

One implication of this scaling form is that the location 
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FIG. 6: Same as Fig.|^ now in 3d with L=20 and line density 
only up to p = 0.5 due to increased computational costs. The 
data is smoothened by taking running averages of each 10 
subsequent data points. 



of the jump of Pperc from to 1 for finite size L saturates 
at some value He in the limit L —f oo, another that the 
jump width decreases with L. 

This appears not to be the case for splay disorder. As 
one can see from Fig. |HIa) He does not saturate for the 
computationally accessible system sizes {L < 120). From 
the inset one can also see that the width of the transition 
does not decrease. 

For a comparison, we use the following model con- 
sisting of random walks embedded into a box of size 
L X L X H. Since the lines typically change their di- 
rection only in the vicinity of the boundaries we consider 
ballistic random walks. They start from a random point 
at the basal plane and evolve towards the top level with 
a random tilt angle taken from the same distribution as 
for the splayed columns. When the random walk meets 
the system boundary it is assigned a new tilt angle. This 
is repeated until the path hits the top of the system. 
This model corresponds to the high density limit where 



FIG. 7: (a) Definition of the winding angle of two fiux lines: 
For each z-coordinate the vector connecting the two lines is 
projected onto that basal plane, z — gives the reference line 
with respect to which the consecutive vectors for increasing 
z-coordinate have an angle (/f>(z). Once (j}{z) > 2n the two 
lines are said to be entangled, (b) Top : A configuration of 
six entangled lines with splay disorder. Bottom: The projec- 
tion of the line configuration on the basal plane, defining a 
connected cluster, (c) Six entangled lines with point disorder. 



all defects are occupied, i.e. B = B^. 

We calculate the mutual winding angles and check for 
entanglement as before. By comparing Fig.lSJa) and Fig. 
IHJb) one observes no qualitative differences in the depen- 
dence of the percolation probability Pperc from the sys- 
tem height H. No saturation is found for percolation 
height He as the system width is increased, though the 
growth is extremely slow. According to our numerical 
results the transition height appears to grow asymptoti- 
cally like He ~ (logL)° with a close to 0.5 (Fig. [TU^ . 

The inset of Fig. |51 shows the relative size of the perco- 
lating bundle Upere/ntot for the data of Fig.|SJa). It tends 
to zero for increasing system size also for H > He ■ This is 
in contrast to what happens for conventional percolation 
where Uperc/ntot — exactly at the threshold and then 
increasing as ~ {H — He)^ with the critical exponent (3. 
Fig. |51 shows the distribution P{k) of the number of lines 
k with which a given line is entangled (i.e. wind around 
it by an angle larger than 2tt). One sees that the proba- 
bility distribution P{k) is almost identical in the cases of 
splay disorder and ballistic random walks, whereas P{k) 
for point disorder is much narrower. 

We conclude that in the case of splayed columnar disor- 
der the spanning bundles are formed mostly by few lines 
which bounce from the boundaries. Thus, there is no 
entanglement percolation transition in the case of splay 
disorder and no giant entangled cluster in the thermody- 
namic limit. 

This behavior changes when the splay disorder is per- 
turbed with attractive or repulsive point defects. Here 
we discuss results for relatively weak attraction. The en- 
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FIG. 8: Percolation probability as a function of the height 
of the system H: (a) splay disorder and (b) ballistic random 
walk. Insets: data shifted with He- For both the transition 
width does not change, there is only a shift of the jump Hc{L). 



ergy costs per unit length for a fraction / of the bonds 
are now set to w = 0.5 whereas the rest of the bonds are 
as before. Fier. 1111 shows the behavior of the percolation 
probabihty corresponding to / = 0.1. It indicates now 
the presence of an entanglement transition following the 
conventional percolation scenario (Eq. {T))). The transi- 
tion height He saturates in the limit L oo (Fig. I1U|) 
and the data obey finite size scaling with the correlation 
length exponent = 4/3 of the 2d percolation problem. 

The fraction / tunes the typical length that a given 
line stays on one defect. When this length exceeds the 
one needed to traverse the system in the lateral direction 



FIG. 9: A plot with logarithmic j/-axis of distribution of the 
degree of lines' entanglement P{k). k is defined as the number 
of lines with which a given line is entangled. The following 
data sets are shown: flux lines in splay disorder with L=100, 
H=115 (ovals); ballistic random walks with L=100, H=39 
(squares); flux lines in point disorder with L=100, 11=48 (tri- 
angles up); flux lines in splay disorder with fragmentation 
L=100, H=70 (triangles down). The system heights are cho- 
sen so that the system is in the vicinity of the entanglement 
transition. Inset: strength of the percolating cluster of fiux 
lines in splay disorder is plotted in a linear scale vs. system 
height H for different system widths, L=30 (ovals) ... 110 
(filled ovals). 
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FIG. 10: A loglog plot of estimated critical system heights 
He vs. logarithm of the system width L. From top: fiux lines 
(FL) with free and fixed (fix) starting points in splay disorder; 
flux lines in splay disorder with fragmentation (fragm.) and 
in point disorder; ballistic random walks (RW), solid line is a 
flt to few last data points. 
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H 



FIG. 11: Percolation probability as a function of the height 
of the system H for splay disorder with fragmentation. Inset: 
corresponding scaling plots, Eq. O. Fragmentation causes 
the sharpening of the transition width with growing L and 
data scales with the 2d percolation correlation length expo- 
nent 1/ — 4/3. 

one observes the properties of pure splay disorder. Con- 
sequently as / is decreased one needs increasingly larger 
system sizes in order to observe the percolation transi- 
tion. The strength and the nature - repulsive or attrac- 
tive - of the perturbing point disorder only changes the 
crossover size of the system as long as the perturbation 
is strong enough to induces lines to change the defects. 
However, due to computational limitations we did not 
attempt to find whether there is a threshold value for 
strength of point disorder below which the lines cannot 
be entangled also in the infinite system size limit. 

VI. SUMMARY AND DISCUSSION 

The ground state of a multi-line system with splay dis- 
order has four different phases as a function of the system 
height H. For small values, there is first a cross-over from 
trivial behavior to the single- line ground state roughening 
as lines start to jump between columnar defects. For in- 
termediate H, one observes the collective regime where C 
depends on dimension but not on the splay energy distri- 
bution. This arises when the lines transverse wandering 
becomes of the same order of magnitude as the average 
line distance. Finally, there is the cross-over to satu- 
ration for L finite. The exponent values are a random 
walk-like C,2d. = 1/2 and a ballistic-like Csd = 1- The 
cross-over between the collective and single-line scalings 
is visible in our 2d data. In 3d, this is not possible for 
numerical restrictions. Note the fact that the 3d single 
line exponent (which varies with the energy distribution) 
is smaller than the collective one. 



We have also considered the separability and entangle- 
ment, to look at the other aspects of these systems as H 
is varied. The former measures the collective aspects of 
the line configuration, which are most pronounced for z 
small. In contrast to the 3d point disorder case, the wind- 
ing of lines around each other is suppressed, leading to 
the absence of entanglement in the thermodynamic limit. 
In the context of our current model, this means that for 
low line densities {B < B^) the enhanced pinning ex- 
pressed by increased critical current, is not a collective 
effect. 

The results for the splay disorder are not expected to 
change qualitatively when the flux lines are allowed to 
enter the homogenous bulk between the defects. As dis- 
cussed in end of SecElin this case there will be a fraction 
of lines staying in the bulk throughout the sample. Since 
these lines have no transverse fluctuations they reduce 
the average roughness by a constant factor. 

Our numerical results suggest that inserting point-like 
defects into splay disorder may even further enhance the 
flux line pinning as the entanglement transition of flux 
lines is recovered. For large line densities {B > B^) - 
where experiments Q imply entanglement - one could 
study the role of additional point disorder. This brings 
new complications, already for single lines. In the case 
of one attractive columnar defect competing with point 
disorder in 2d one finds only a localized line, whereas 
in 3d there a localization-delocalization transition |2lj |. 
Another study has found for a mixture of many columnar 
defects and point disorder sub-ballistic behavior p^. 

Lidmar et al. found by mixing in point disorder that 
in 2d splay disorder always dominates on the long length 
scales whereas in 3d it seems that strong enough frag- 
mentation leads to point disorder behavior ^3 • Here we 
briefly present a few possible scenarios associated with 
tuning the strength of point disorder with fixed strength 
of splay disorder: (i) At large enough system sizes point 
disorder will always dominate, (ii) There is a crossover 
from pure splay disorder behavior (no entanglement) to 
point disorder behavior, (iii) There is a third regime 
where lines take an advantage of both splay and point 
disorder leading possibly to more efficient entanglement 
than in the case of pure point disorder. This would be 
due to the large displacements along splayed defects. One 
can expect that the response of the system to point disor- 
der perturbations depends on the line density. It would 
be interesting to find the full phase diagram including 
the optimal parameters from the point of view of entan- 
glement. 
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